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In a previous Letter [Phys. Rev. Lett. 96, 123401 (2006)] we have shown by means of three- 
dimensional particle-in-cell simulations and a simple rigid-sphere model that nonlinear resonance 
absorption is the dominant collisionless absorption mechanism in the intense, short-pulse laser cluster 
interaction. In this paper we present a more detailed account of the matter. In particular we 
show that the absorption efficiency is almost independent of the laser polarization. In the rigid- 
sphere model, the absorbed energy increases by many orders of magnitude at a certain threshold 
laser intensity. The particle-in-cell results display maximum fractional absorption around the same 
intensity. We calculate the threshold intensity and show that it is underestimated by the common 
over-barrier ionization estimate. 



PACS numbers: 36.40.Gk, 52.25.0s, 52.50.Jm 



I. INTRODUCTION 

Intense laser-matter interaction provides a route to 
create energetic particles (e.g., electrons, ions, and pho- 
tons) using table-top equipment. Clusters, possessing the 
transparency of gas targets and the high charge density of 
solid targets, proved to be very efficient absorbers of laser 
light. Their small size, compared to laser wavelength and 
skin depth, avoids reflection of the laser beam at the clus- 
ter surface as well as the loss of hot electrons into the cold 
bulk. In fact, almost 100% absorption of the laser light 
was reported in experiments with rare _ga_s_ clusters [l| . 



Useful reviews on the subject are Refs. liia. 

Upon irradiation of the rare gas clusters by intense 
laser light, electrons first absorb energy and leave their 
"parent" ions. This is known as inner ionization, meaning 
that the electrons are still bound to the cluster but not 
necessarily to their "parent" ions. The total electric field 
(i.e., laser plus space charge field) inside the cluster leads 
to inner ionization up to high charge states not possible 
with the laser field alone (ionization ignition As 
the laser intensity during the pulse increases, these elec- 
trons absorb energy from the laser field and may leave 
the cluster leading to the positive charging of the cluster 
known as outer ionization. Thus outer ionization leads to 
a non-neutral plasma. With the increasing outer ioniza- 
tion, the restoring force of the ions counteracts ionization 
ignition so that the latter stops at some point. The net 
positive charge left behind finally explodes due to the 
Coulomb repulsion and hydrodynamic pressure, leading 
to the conversion of electron energy into ion energy. Typ- 
ically MeV ions and keV electrons |3i S 13 are measured 
in experiments. 

It is clear from the described scenario that the un- 
derstanding of the relevant mechanism(s) of laser energy 
absorption leading to the heating of cluster electrons and 
outer ionization is of great importance for the develop- 
ment of a complete theoretical description. 

Laser energy absorption by electrons proceeds either 
through resonances (linear or nonlinear) or through non- 
adiabaticities (all possible types of collisions). All of 



these processes lead to dephasing of the current density 
with respect to the laser field, which, according to Poynt- 
ing's theorem, is a prerequisite for absorption. CoUisional 
absorption via collisions of electrons with ions are of mi- 
nor importance at near infrared wavelengths ~ 800 nm 
or greater |^ 0, ^| whereas it is the dominant 
absorption mechanism at short wavelengths [TsL Q, , 
not studied in this paper. The finite size of the clusters 
suggests that "collisions with the cluster boundary" may 
be responsible for the energy absorption. However, this 
viewpoint is misleading, as will be shown in Sec. IIIBI 

During the expansion of the ionic core, the decreas- 
ing charge density p{t) leads to the decrease of the Mie 
plasma frequency, WMio(i) = ■\/47rp(t)/3 (atomic units 
h — m — —e = 47reo = 1 are used unless noted other- 
wise). For very short near infrared laser pulses ujuieit) 
cannot meet the linear resonance 
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unless the cluster has sufficiently expanded (typically af- 
ter a few hundred femtoseconds). Linear resonance 
well understood in theory, experiments, and simulations 
[lllIlEEHEl|2i|23|,|23, is thus ruled out for 
very short pulses or during the early cycles of a long pulse 
laser-cluster interaction where ion motion is negligible. 
In this case, nonlinear resonance (NLR), whose origin 
lies in the anharmonicity of the cluster potential, turns 
out to be the dominant collisionless absorption mecha- 
nism. In fact, for very short linearly polarized (LP) laser 
pulses, it was clearly shown 25] that essentially all elec- 
trons that contribute to outer ionization pass through the 
NLR, which was unequivocally identified as the collision- 
less absorption mechanism in the absence of linear res- 
onance. The eigenfrequency uj[r(t)] of a (laser-) driven 
oscillator in an anharmonic potential, being dependent 
on its excursion amplitude r{t) (or the energy), may dy- 
namically meet the NLR 
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Due to the many body nature of the interaction, the 
identification, the separation, and the interpretation of 
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the absorption mechanisms in molecular dynamics or 
particle-in-cell (PIC) simulations are often difScult. Re- 
cently, a method of identification of the NLR in many- 
body simulations of rare gas clusters has been proposed 
|25| . The possible importance of NLR was also mentioned 
or discussed previo usly l2^ l27i l2^ l2^ . The rigid sphere 
model (RSM) pi IM Isif where electrons and ions are 
modeled by homogeneously charged rigid spheres oscillat- 
ing against each other is clearly an oversimplification of 
a real many-particle system such as a cluster. However, 
it proves useful for estimating the order of magnitude of 
the absorbed energy as well as for the calculation of the 
laser intensity where energy absorption is most efhcient, 
as will be shown in the present work. Moreover, it pro- 
vides physical understanding and clearly displays NLR 

Mm 

The heating of cluster electrons in circularly polarized 
laser fields has not yet received much attention, at least 
theoretically. Experiments with rare- gas clusters show al- 
most no effect on the x-ray emission [S^, 133, ^'^^ 
energy distribution when laser light of different ellip- 
ticity is used. Theoretically, circular polarization is par- 
ticularly interesting because the above mentioned "colli- 
sions with the cluster boundary" are strongly suppressed 
in this case. Hence one may expect energy absorption 
being less efficient. NLR, on the other hand, occurs in 
both cases, and, in fact, the energy absorption turns out 
to be equally efficient. Moreover, the absorbed energy 
compares well with the RSM predictions. 

The outline of the present paper is as follows: in 
Sec. HTKl we briefly review the NLR and the RSM. In 
Sec. Ill Bl the RSM is extended to circular polarization, 
where NLR is observed as well. In Sec. IhTI the RSM 
threshold intensities are calculated. Section llTTl is devoted 
to the PIC results for both linear and circular polar- 
ization. Finally, we summarize our results in Sec. IIVI 

Throughout this paper we use n = 8-cycle laser pulses 
of near infrared wavelength A = 1056 nm and a fixed clus- 
ter radius R — 3.2 nm unless stated otherwise. NLR is 
a robust phenomenon that — qualitatively — is insensitive 
to cluster and laser parameters. 



at the center of the ion sphere. In an oscillating laser 
field, the homogeneously charged electron sphere oscil- 
lates back and forth against the positively charged ion 
sphere. 



A. NLR in a linearly polarized laser field 

The equation of motion of the electron center of mass 
in a LP laser field, polarized along x, can be written as 



dV f , , Ei{t) 



(3) 



where f = x/R is the excursion of the electron sphere, 
normalized to the cluster radius R, r = \f\, t — ujit is 
the normalized time in units of the laser period, and the 
dimensionless electrostatic restoring force g{r) is 



r-'4 + T2 0<-<2 (4) 
4 r > 2 ' ' 



The first term in the upper line of is the linear force 
when the displacement of the electron sphere is small, 
the next two terms are the nonlinear terms which ap- 
pear due to the partial overlap of the electron cloud 
with the ion cloud. The term in the lower line of 
Q is the Coulomb force between the separated electron 
sphere and ion sphere. The quantity on the right hand 
side of © is the normalized driver strength, which is the 
quiver amplitude E\{T)/ujf of a free electron in the laser 
field ^i(t) normalized to the cluster radius R. Earlier 
work [23 showed that absorption of laser energy in the 
RSM is characterized by a threshold driver strength be- 
low which absorption is negligible (harmonic regime) and 
above which absorption is almost constant. Figure ^ 
shows this threshold behavior for (tJMie/wi)^ = 10/3 and 
a, n = 8-cycle sin^-pulse Ei{t) — Eo sin^(T/2n) cos(r) for 
< r < 2mr. The dashed line is the absorption corre- 
sponding to a driven, purely harmonic oscillator 
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II. NONLINEAR RESONANCE IN THE RIGID 
SPHERE MODEL (RSM) 



in the laser field Ei{t). The laser energy absorbed by a 
single electron in a n-cycle laser pulse of period T is 



The ion motion can be neglected in the study of energy 
absorption in very short laser pulses. Thus the ions just 
form a static, positively charged background of spheri- 
cal shape. For not too high laser intensity, the collective 
motion of the electrons can also be approximated by a 
homogeneous, rigid sphere of negative charge. In the 
simplest case the radii of ion and electron sphere are as- 
sumed to be equal. In a more realistic model the electron 
cloud expands |23|. However, the method of identifica- 
tion of the NLR used in this paper is independent of the 
degree of electronic expansion. The center of mass of the 
electron-ion system is, in good approximation, located 
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v{t) ■ Ei{t)dt. 
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Solving |SJ analytically for the velocity v{t) and integrat- 
ing ® one finds for the absorbed energy per electron 
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1 - cos (wMie^T)] (7) 



where fw^;, 



(l-|-l/n)VF[<,-(l-l/n)V 



('^Mic~^i ) (f°^ sufhciently large n) was used. The ana- 
lytical estimate {TJ is plotted in Fig. together with the 
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FIG. 1: (color online). Laser energy absorption vs driver 
amplitude in the rigid sphere-model for (iOMie/i^i)^ = 10/3, 
a n = 8 cycle laser field Ei{t) = i5o sin^(r/2n) cos(r), and 
cluster radius R = 3.2 nm. Within a narrow field strength 
interval (here ~ 0.5-1) the absorbed energy per particle (solid 
line) increases by many orders of magnitude. The dashed 
line represents the absorbed energy Q by a purely harmonic 
oscillator driven by the same laser field. 



absorbed energy obtained from the numerical solution of 
One sees that the absorbed energy jumps by many 
orders of magnitude to a higher value after crossing a 
threshold driver strength. The higher the cluster charge 
density is, the higher is this threshold driver strength. 
The rigid sphere model shows this behavior of efficient 
absorption above the threshold driver strength at all clus- 
ter charge densities independent of the linear resonance 
condition p = Spc, contrary to the nanoplasma model 
[iGj . Since the rigid sphere model does not necessarily re- 
quires expansion of the cluster for the efficient absorption 
of laser energy, it permits us to understand the behavior 
of energy absorption and the underlying mechanism for 
very short laser pulses. 

Equation can be formally rewritten as 



d r 



Riuf ■ 



(8) 



Equation ^ yields the instantaneous, scaled effective fre- 
quency squared 



Ei{t 



r{T) 



r{T) 



r(r) 



(9) 



which passes through unity at the NLR (0). The right 
hand side of © is the restoring force divided by the ex- 
cursion of the electronic cloud, which in the case of har- 
monic motion would be the square of the characteristic 
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FIG. 2: (color online). Typical behavior of {LOesir) /uiif (red, 
upper solid line) vs laser cycles above the threshold driver 
strength for (ojMic/i^i)^ ~ 40/3. Here we take Eo/Rcof ~ 7.48 
corresponding to a laser intensity ~ 2.5 x lO^^W/cm^, a n = 
8-cycle sin^-pulse -Ei(r) = Eo sin^(r/2n) cos(r) of wavelength 
A = 1056 nm, and a cluster radius R — 3.2 nm. Excursion x/R 
(green, middle solid line) and energy of the electron sphere 
Etot/R^(^'\ (blue, lower solid line) are included in the plot. 
Outer ionization (i.e., Etot/R^iOi > 0) and occurrence of NLR 
[LOca{T) /loi^ = 1 always coincide (dashed vertical line). 



frequency. Figure |21 shows a typical example of the tem- 
poral behavior of [wcff (T)/aJi]^ above the threshold driver 
strength for {uoyiic/^^xY — 40/3. Since (wmic/^^i)^ — 
13.32, [LJcfi{'T)/ijJ\^ starts at this value and drops with in- 
creasing driver strength. It passes through unity at the 
time indicated by the vertical line, and it is exactly at 
that time where the electron sphere becomes free (outer 
ionization) . This incidence is clearly visible from the en- 
ergy of the electron sphere, which passes through zero 
and the excursion as well. Outer ionization and occur- 
rence of NLR happens for all driver strengths above the 
threshold whereas the resonance is never met below the 
threshold. Since the amplitude of the excursion of the 
electronic sphere depends upon the driver strength, the 
excursion amplitude should also be large enough so that 
the NLR is passed. The decrease in the effective fre- 
quency with the increase of the amplitude of excursion 
of the electronic sphere in the force field © can be un- 
derstood by analyzing its motion in the corresponding 
anharmonic potential 



V{r) = uJMijR' X 



— — n<r<2 

2 16 160 u r ^ z 

f-- r>2. 

or — 



(10) 



The period T of oscillation of the electronic sphere in the 
potential V{r) can be approximated by a perturbation 
series [sjl as 



(-1)'= 9* 



k\ 



be'' 



[5V{r)\' dr 



(11) 



Here, i?tot is the total energy of the electronic sphere 
in the ionic field, Vb(r) = a'^jj,i?^r^/2 is the harmonic 
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FIG. 3: (color online). Effective frequency ujc«/i^i vs the 
excursion amplitude f of the electronic sphere for different 
cluster charge densities p/pc ~ 10-40. The solid lines are 
computed from the numerical solution of The dashed 

lines are the corresponding analytical approximations from 
IIH for k up to 6, and the dashed-dotted lines are from 1121 . 
For the charge density p/pc = 40 one expects NLR (horizontal 
dashed line) to occur at the excursion f ~ 2.88. 



oscillator potential and dV{r) = ai?,'^r'^/3 + (3R^r^ /5 is 
the perturbation to the potential with a = — 9w^j^/16i? 

i^/32R^. The effective frequency is then 



and (3 

Woff — 2tt/T. For the excursion r < 2 and cluster ra- 
dius R = 3.2nni (~ 60.4a.u.) we can consider \P/a\ ^ 1 
and the approximate potential SV{r) ~ aR^r^ /3. Cor- 
rections up to fc = 6 yield T = Tq + Ti -\ Tq. The 

unperturbed period is Tq = 27r/cjMic and the succes- 
sive corrections are Ti = ci/{3uj^-J, T2 = 02/(3^^;^), 
Ta = C3/(9<J, r4 = c,/{72loI!J, T5 = c,/i5Au;l^J, 
Tq = C6/(36a;ifiJ with ci ~ -8a(2£;tot)l/^ C2 ~ 
5TTa^Etot, C3 ~ -28A5a^{2Etot)^/^, C4 ~ 385^4^^^^, 
C5 ~ -318.6a5(2£;tot)^/^ and cg ~ 97.89 x 2^a^E^^^, re- 
spectively. 

Neglecting the higher order term of g{r) for r < 2, a 
simpler approximation to the effective frequency can be 
derived from by replacing the excursion r with the 
excursion amplitude r: 

uj,s{f) ~ t^Mic (1 - 9f/16)'/' ~ WMic (1 - 9f/32) . (12) 

Figure O shows the effective frequency vs the excursion 
amplitude f of the electronic sphere in the potential H1(J|) 
for various cluster charge densities p/pc — 10-40 as calcu- 
lated from the numerical solution of together with the 
approximations (|ll|l and H12() . The effective frequency as 
calculated from (|11|) shows good agreement below the ex- 
cursion f < 1.5 and low charge densities (e.g., p/pc < 10). 
For higher charge densities more corrections [large num- 
ber of k values in Hll|l ] are needed, which are very much 
cumbersome to calculate. Although (|12|l fits well with the 
numerical solution, lOcS becomes negative when f > 32/9. 
However, the results using H12|l agree with the exact ones 



in Fig. O up to p/pc = 30 for the excursions of interest, 
i.e., up to the point where the electronic sphere under- 
goes NLR. The variation of frequency with the excursion 
amplitude as shown in Fig. |31 explains why a threshold 
driver strength is required for an appreciable laser energy 
absorption as well as for the crossing of the NLR in the 
Fig. 121 only a driver exceeding a certain threshold field 
strength will lead to excursions compatible with the NLR 
condition lUcs — lu\. 



B. NLR in a circularly polarized laser field 

Clusters in a circularly polarized (CP) laser field re- 
ceived less attention in the literature. It is not known 
a priori how the outer ionization and energy absorption 
by clusters depend on the laser polarization. In laser- 
atom interaction the laser polarization has dramatic ef- 
fects: since in CP the free electrons do not return to their 
parent atom all the atomic effects relying on rescatter- 
ing such as high-order harmonic generation, high-order 
above-threshold ionization, and nonsequential ionization 
are strongly suppressed. In the context of clusters, the 
study of the absorption efficiency as a function of the 
laser polarization can help to discriminate among differ- 
ent absorption mechanisms. For instance, if laser en- 
ergy absorption was indeed due to "collisions with the 
cluster boundary" it would be suppressed in CP because 
the electrons mainly swirl around parallel to the "clus- 
ter boundary" rather than crossing (and hence colliding) 
with it. However, as we will show, the absorption of laser 
energy is largely independent of the laser polarization, 
thus ruling out "collisions with the cluster boundary" as 
a meaningful absorption mechanism. 



1. Two-dimensional rigid sphere model 

Let us first extend the RSM to CP. In a CP laser field 
with electric field components in x- and y-direction, the 
equation of motion for the electronic sphere in the rigid 
sphere approximation of a cluster can be written as 



r 



1 



Efir) 



(13) 



Here, r^ = x/R, ry = y/R, 



ry , and 



Efir) 



Eo_ 
V2 



sin^(T/2n) cos(t) 



V2 



sin (T/2n) sin(T). 



(14) 
(15) 



One identifies E^ = g{r)rx/r and Ey — g{r)ry/r as the 
two components of the restoring force in (|13|l . Note that 
we have divided the electric field components by a factor 
V2 so that the ponderomotive energy Up = Eg/Auif, (i.e., 
the time-averaged quiver energy of a free electron in the 
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Laser cycles 

FIG. 4: (color online). Typical behavior of (aJeft(T)/o;i)^ 
(green, top left solid line) vs laser cycles above the thresh- 
old intensity ~ 3 x lO^^W/cm'^ for circular polarization, 
(i^Mic/i^i)^ = 40/3, and cluster radius R = 3.2 nm. Here, 



i.2, corresponding to a laser intensity ~ 5 x 
-cycle sin^-pulse with components Ei{t) 



Eo/Rujf ~ 
lO^'^W/cm^ 

£;osin2(r/2n)cos(r)/^/2, E^{t) = E^o sin^(T/2n) sin(r)/v^, 
and wavelength A = 1056 nm. Excursion r (blue, middle 
left solid line) and energy Etot/R^i^i (red, bottom left solid 
Une) of the electron sphere are also plotted. The energy 
Etot/R'^^i is scaled down by a factor 10 to display within 
the excursion and the frequency range. Outer ionization (i.e., 
Etot/R^(^i > 0) and occurrence of NLR [LJcff{r)/ujif — 1 
always coincide (dashed vertical). 



laser field) is the same as in the LP case with the same 
£^0 (otherwise Up would be a factor of two higher in the 
CP case). The square of the effective, time-dependent 
oscillator frequency in the CP laser field can be written 
as 



^cff(T) 



'^y-^x 



r(r) 



(16) 



which has the same form as in the LP case 0. Earlier, 
in Fig. [3 it was shown that NLR and outer ionization 
in the RSM only occur when a threshold laser intensity 
is crossed. The same is true for the occurrence of NLR 
with CP light. Figure 01 shows the temporal behavior 
of [uJes{T)/oJi]'^ above the threshold driver strength for a 
n = 8-cycle CP laser pulse of wavelength A = 1056 nm. 
Here the cluster charge density is 40 times over-critical, 
i.e., (wMic/wi)^ = 13.32 at which (woff (t)/wi)^ starts and 
drops with increasing driver field during the pulse. The 
NLR [ll>cs{t)/uji]'^ = 1 is passed at the time indicated by 
the vertical line. As in the LP case in Fig. 12 the elec- 
tron sphere is set free at the time the NLR is passed: the 
energy of the electron sphere passes through zero, and 
the excursion sharply increases to a high value. It is also 
clear from Fig. ^ that once the electronic sphere is set 
free, the frequency drops to zero, and the total absorbed 
energy remains positive. A zero effective frequency im- 
plies an infinite period, i.e, the electron sphere does not 
return to the ion sphere. The main difference between 
Fig. ^ and Fig. [21 is that in the case of the CP laser field 



the decrease of the effective frequency is smooth (i.e., 
no oscillations) since the electric field vector rotates but 
its absolute value remains constant. As a consequence, 
the electron sphere spirals out, staying away from the 
potential center where lOcS = iouie- For LP instead, the 
electron sphere is driven through the origin and hence the 
effective frequency undergoes oscillations before it drops 
to the resonance value, as visible in Fig. 13 NLR is clearly 
identified in both cases. 



C. Prediction of threshold intensity for the NLR 

NLR occurs above a threshold driver strength. Be- 
yond this driver strength the rigid electron sphere gains 
laser energy which is many order of magnitude higher 
than below the threshold (see Fig. In an open po- 
tential such as (|10|l the electron sphere is detached from 
the ion sphere above the threshold driver strength, i.e., 
appreciable energy absorption and outer ionization occur 
simultaneously. We would like to remark that in closed 
potentials NLR occurs as well, as discussed in Ref . 28] . 

The threshold driver strength can be estimated. The 
dimensionless potential U{r) = l^(r)/(w^;(,i?^) of the 
electron sphere in the ionic field can be written as 



H 3r'^ 

U(r) = \ , 

^ ' 2 16 160 



r < 2. 



(17) 



Application of a static electric field Eq (corresponding 
to the peak field strength of a low- frequency laser field), 
suppresses the potential in one direction by the amount 
REq r. The effective potential seen by the electron sphere 
is i7off(r) =U{r)-Eor with = En / {tol^^.R) . The po- 
tential barrier vanishes if U' {rv)—Eo = and U"{rv) — 0, 
leading to Tv = 1, and the NLR threshold intensity is es- 
timated to be 



rVBA 
^th 



5 
32 



(18) 



We call this the vanishing-barrier-approximation (VBA) . 

In atomic ionization the so-called over-the-barrier ap- 
proximation (OBA) or Bethe-rule js^ allows to esti- 
mate at which electric field strength a certain atomic 
charge state dominates. If one applies the OBA to the 
RSM one obtains the two equations Ues{rh) = and 
U'{rt,) — i?o = with rb the barrier location. One finds 
Tb ~ 1.613. This gives the OBA threshold intensity of 
the NLR 



OBA 

th 



(19) 



which underestimates the numerically determined 
threshold intensity, as we will show now. 

Figure O shows the threshold intensity as a function of 
the cluster charge density p/ Pc = 10-160 (corresponding 
to average charge states ~ 1 — 16 for Xenon). The numer- 
ically determined threshold intensities for LP (j^l and CP 
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FIG. 5: NLR threshold intensities in the RSM vs the cluster 
charge density p/ Pc for linearly (LP) and circularly (CP) po- 
larized laser fields. Results from the full numerical solution of 
^ with LP (lower solid), 1131 CP (top solid), the vanishing 
barrier approximation (VBA-LP, middle solid) 11811 . and the 
VBA corrected for CP (upper dashed) are shown (see text 
for a discussion) . The over-the-barrier approximation (OB A) 
1191 1 (lower dashed) underestimates the exact threshold inten- 
sities. 

(|13|l laser light show that when the cluster charge density 
is low, the NLR occurs almost at the same value of the 
threshold intensity, irrespective of the polarization. As 
the charge density increases the NLR threshold intensity 
appears to be higher for CP than for LP. The VBA l(TS|) 
of the threshold intensity is in good agreement with the 
numerical result for LP whereas the OBA H19|l under- 
estimates it. This fact might be related to the recently 
observed "enhanced saturation intensities" in the ioniza- 
tion of finite size systems such as Cgo (see, e.g., l40l | 
and references therein), indicating that the latter might 
neither be a many-electron nor a quantum effect but just 
due to the finite size of the target. 

The difference of the threshold intensities for LP and 
CP is due to the definition of the CP field (fn|l , ifTCjl where 
a factor 2~^/^ has been introduced in order to render the 
ponderomotive potential equal for LP and CP. However, 
for the threshold intensity it is the electric field (or the 
intensity) that matters, not Up. For a given Eq the laser 
intensity is Iq = Eq in the LP case but only Iq/2 for CP. 
Therefore, the upper black, dashed line in Fig. |S1 shows 
the VBA threshold intensity multiplied by a factor of two, 
which is in good agreement with the numerical results for 
the CP laser field at higher charge densities. 

So far we have studied the NLR absorption of laser 
energy in a simplified model system assuming an anhar- 
monic potential generated by the ions in which the ho- 
mogeneous and rigid electron cloud moves. In reality, 
the potential builds up during the interaction with the 
laser pulse because of ionization. The delicate interplay 
of inner ionization, energy absorption by various mecha- 
nisms, and outer ionization can be simulated using meth- 
ods such as PIC or molecular dynamics. Previous work 
[25| . studying LP short laser pulses, showed that NLR 



FIG. 6: (color online). Kinetic energy (KE) and total ab- 
sorbed (TE) energy per electron in units of Up vs laser in- 
tensity for charge density 40 times the critical density and a 
Xeieii cluster (N = 5498 electrons, mean ionic charge state 
~ 3.4, and a cluster radius R = 3.2 nm). Solid lines are 
PIC results, dashed lines are the RSM results for a n = 8 
cycle pulse Ei{t) = Eosm'^{cjit/2n) cos{LJit) of wavelength 
A = 1056 nm. 

(0) can be clearly identified in such simulations as well. 



III. NONLINEAR RESONANCE: 
PARTICLE-IN-CELL RESULTS 

In this Section we present results obtained from three- 
dimensional PIC simulations. The cluster is exposed to 
n — 8-cycle sin^-pulses Ei{t) — Eq siii^ {uj\t / 2n) cos{uJit) 
of near infrared wavelength A = 1056 nm, i.e., the total 
pulse duration is 28 fs. Since motion of ions does not 
play an important role during the entire pulse the ions 
are assumed fixed, which ensures a well defined, constant 
Mie frequency wmic- A PIC electron has the same charge 
to mass ratio as a "real" electron, that is, e/m = — 1 
in atomic units. Each PIC electron moves under the in- 
fluence of the external laser field and the space charge 
field Esc = — V<I>(r,t) due to the potential <i>(r,t) that 
is created by all charges (mapped to the numerical grid) . 
Hence the equation of motion of the i-th PIC electron is 

r, + Esc{n,t) = -Ei{t). (20) 

Equation (|2U|) is solved self-consistently for all PIC elec- 
trons. Clearly, Esc{ri,t) depends on the position of all 
other particles ^ i as well. The PIC simulation starts 
with the neutral cluster configuration so that Esdvi, 0) = 
0. 



A. Results for linear polarization 

FigureElshows PIC results for the kinetic energy (KE) 
and the total energy (TE) absorbed per electron in units 
of Up, vs the peak laser intensity. One sees that the 
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absorbed energy per electron is on the order of Up. How- 
ever, the absorbed energy is nonhnear in Up and dis- 
plays a maximum before it drops because of the satu- 
ration of outer ionization. The depletion is due to the 
fact that at a given intensity most (if not all) electrons 
are removed from the cluster (complete outer ionization) . 
Further inner ionization would be required to generate 
"fresh" electrons that could continue to absorb energy. 
The maxima in the PIC absorption curves are located 
close to the threshold intensity predicted by the RSM. 
The total absorbed energy around the maximum of the 
PIC absorption curves (RSM as well) is on the order of 
2Up = 5-6 keV which has been also reported in experi- 
ments of intense laser clusters interactions . 

In PIC simulations no sharp intensity threshold exists 
since each PIC electron sees its own time-dependent field 
(space charge field plus the laser field). Therefore sharp 
jumps (as seen in the RSM, e.g., Fig.P) are absent. 

Figure |7| compares the kinetic and potential energies 
of a single PIC electron and the electronic sphere in the 
RSM. The energies are plotted vs the excursion x/R 
(along the direction of the laser field). The peak laser 
intensity is Iq — 2.5 x 10^^ W/cm^. We have checked 
that there are many PIC electrons that leave the cluster 
similarly to the one presented in Fig. [7| The potential 
and the kinetic energy gained by the electron sphere in 
the RSM is divided by the total number of PIC elec- 
trons so that a quantitative comparison is possible. The 
main difference between PIC and RSM results is that in 
the PIC simulation all the electrons see initially a zero 
potential since the potential builds up from the neutral 
cluster configuration in the course of the interaction with 
the laser field while in the RSM the electronic sphere os- 
cillates in a prescribed potential. Those PIC electrons 
that stay long inside the cluster experience almost the 
full ionic potential which then is similar to the RSM po- 
tential. 

As mentioned earlier, the absorption of energy by a 
PIC electron depends upon the self-consistent potential 
which develops during the laser pulse due to outer ion- 
ization. As a result different PIC electrons move along 
different trajectories, "see" a different potential, and thus 
are set free at different times. Figure |S1 shows a snapshot 
of the collective potential ^{r,t) at time u)\t/2n — 4.81 
for a cut at z = (and various y throughout the clus- 
ter). The most lower curve represents the potential for 
y = z = Q. The red circles represent the total energy 
Etot,i = r'f{t)/2 — $(ri,t) of individual PIC electrons 
located within the simulation box at that time. 

One observes that many PIC electrons are accumulated 
near the bottom-left of the potential well at this time. 
These PIC electrons with i?tot,i < remain bound since 
the laser field amplitude is already decreasing from the 
4th cycle onward. The kinetic and the total energy of 
a PIC electron which leaves the potential well when its 
excursion becomes x/R ~ 2.88 at the same time is also 
shown. This excursion approximately satisfies the NLR 
condition ujas/i^\ = 1 for the charge density p/ Pc — 40, 
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FIG. 7: (color online). Comparison between the kinetic and 
the potential energies of a PIC electron and the equivalent 
RSM electronic particle with same charge and mass of a PIC 
electron vs the excursion x/R in the laser polarization direc- 
tion. The escaping PIC electron resembles the equivalent elec- 
tronic particle in the RSM. Peak laser intensity 7o — 2.5 x 10^^ 
W/cm'^, charge density is p/pc = 40. Other parameters as in 
Fig.ini 
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FIG. 8: (color online). Energy and self-consistent potential 
$ from the PIC simulation vs the excursion x/R in laser po- 
larization direction. Yellow (solid gray) lines represent cuts 
of the potential $ at time uji\t/2-K = 4.81 for different y and 
2 = 0. The circles (red) represent the total energy Etax,i of 
individual PIC electrons located within the simulation box 
at that time. The dashed lines represent the kinetic (green, 
dashed gray) energy (KE) and potential (dashed black) en- 
ergy (PE) of a PIC electron that is outer ionized when the 
excursion a;/i? ~ 2.88 meets the condition tJcft = ^\ at the 
same time (see Fig.|31for the charge density p/ Pc = 40). The 
peak laser intensity is Jo — 2.5 x lO^'' W/cm^. Other param- 
eters as in Fig. [3 

as identified in Fig. O with the RSM analysis. 

The NLR behavior exhibited by the PIC electron in 
Fig. |S| is not accidental. For the sake of an unequivocal 
and explicit identification of the NLR we now analyze the 
motion of all individual PIC electrons in the same way as 
it has been done with the motion of the electron sphere in 
the RSM in Sec^ Recalling (j^Ol), the equation for the 
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FIG. 9: Effective frequency squared (a)eft,i/wi) , excursion 
Ti/R, and tlie totai energy Etot,i = r?(t)/2 — $(ri,t) for the 
PIC electron of Fig.|7|vs time in laser cycles. The total energy 
becomes positive only when the NLR is crossed (indicated 
by the vertical dashed line). This result resembles the RSM 
result in Fig. |21 when NLR is met. The charge density is 
p/pc = 40, the peak laser intensity is 7o = 2.5 x 10^^ W/cm^. 
Other parameters as in Fig.|3 



effective, time-dependent oscillator frequency analogous 
to H16|l for the z-th PIC electron reads 



2 



[E^{t) + r,{t)] ■ r,{t) E,,{r„t)-r,{t) 



(21) 

Earlier we have mentioned that Esc{ri,t) depends on 
the position of all other particles 7^ i and the simula- 
tion starts with the charge neutral cluster configuration 
i.e., Esc{ri,0) = 0. Hence, a PIC electron "sees" ini- 
tially an effective frequency aJoff,i(0) = 0. The laser field 
disturbs the charge equilibrium and Lu'^g ^ (t) becomes dif- 
ferent from zero, uj'^g j(t) may be even negative in regions 
of accumulated electron density (repulsive potential) . As 
the cluster charges up, (cJeff/'-^i)^ quickly increases be- 
yond unity (where the RSM starts in the first place). 
The start from LUcs^i{0) — 0, the possibility of negative 
^cS i (^) ' ^'^^ three-dimensionality are the main differ- 
ences to the RSM analysis above. Figure El shows the ef- 
fective frequency squared, the total energy Etot,i{i), and 
the excursion Vi/R vs time for the PIC electron whose 
energy history is shown in Fig. [7| We define the time 
when, for a particular electron, Etot,i becomes > as the 
ionization time of that electron. Data are plotted shortly 
before the emission of this particle from the cluster po- 
tential. It is clearly visible in Fig.^that the PIC electron 
is escaped only when the resonance line (lUcs/^i)'^ = 1 is 
passed. Figure|51can be well compared with Fig.[21show- 
ing ionization of the RSM via NLR. 

In the case of PIC simulations the fulfillment of the 
nonlinear resonance condition (woff/t^i)^ — 1 is neces- 
sary but not sufficient for ionization. As the potential 
builds up, the PIC electrons transiently meet the NLR 
condition, and, in fact, some electrons leave the cluster 
at that early stage when the potential is still shallow and 
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FIG. 10: (color online). Snapshots of PIC electrons in the 
frequency vs energy-plane at times (a) t = 1.5, (b) t = 2.0, (c) 
t = 2.5, (d) t = 3.0, (e) t = 3.5, (f) t = 4.0, (g) t = 4.5, and (h) 
t = 5.0 laser cycles for LP, laser intensity 2.5 x 10^^ Wcm~^, 
and {ujMic/(^i)'^ ~ 40/3. Other parameters as in Fig.m The 
radial positions (in units of R) are color-coded. Electrons 
become free upon crossing the NLR, i.e., {uj'^g / ujf , Etot /Up) = 
(1,0). 



the laser field is relatively weak. However, as the poten- 
tial deepens, PIC electrons "dropping" below the energy 
necessary for NLR to occur, behave from then on simi- 
lar to the RSM and may finally escape only by climbing 
up in the potential and hitting the NLR (wcfi/wi)^ = 1. 
During the 0.4 laser cycles plotted in Fig.^lthe {uj^s /^jJi)'^- 
curve displays artificial, short-time scale fluctuations in- 
herent in PIC simulations [s^. However, we checked that 
macroscopic observables such as the absorbed energy or 
the degree of outer ionization are well converged. 

By following the dynamics of the electrons in the ef- 
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fective frequency vs energy-plane we identify the main 
pathway to outer ionization and efficient absorption. Fig- 
ure Elt-h shows the scaled effective frequencies squared 
{lOcs/uj\Y of the individual PIC electrons vs their ener- 
gies Etotit) = rf{t)/2 — <^{ri,t) they would have if the 
driver is switched off instantaneously at t — 1.5, 2, 2.5, 
3, 3.5, 4, 4.5, and 5 laser cycles, respectively. At the 
time when, for a particular electron, Etot becomes > 
the ionization occurs for that electron. The laser inten- 
sity is 2.5 X lO^^Wcm"^, and the pre-ionized cluster is 
40 times over-critical so that [ujyim/^if = 40/3. As is 
clearly visible in Fig. ^| each electron reaches positive 
energy close to the point (w^g/wf , E'tot/C^p) = (1, 0). The 
radial position of each electron is color-coded, indicating 
that outer ionization occurs at radii around 2R. During 
the early time of the laser pulse (Fig. llOb-b) when many 
electrons are still inside the cluster, (weff/t^^i)^ spreads 
over a wide range, starting from the maximum value 
(wMie/wi)^ down to negative values due to the repulsive 
force exerted by the compressed electronic cloud. Note 
that negative values in effective frequency occur mainly 
at early times where most of the electrons are still inside 
the cluster. Electrons with positive but very small i^tot 
and ~ represent low energetic electrons removed 
earlier during the pulse (see Fig.Hnh.b). The occurrence 
of NLR is less clear for these early leaving electrons. As 
mentioned above, these electrons move in a shallow ef- 
fective potential with {uJcs/oJ\f' < 1 when they leave the 
cluster with ease and with rather low kinetic energy be- 
cause the laser intensity is still low at the time of their 
emission. Figures llOb-f show that most of the electrons 
escape from the cluster by passing through the channel 
(uj'^g/ujfjEtot/Up) = (IjO) at radii around 2R. It is also 
visible that more and more electrons are driven to pos- 
itive frequency before they leave the cluster by passing 
through (uj'^g/ujfjEtot/Up) = (1,0). This is so because 
as more and more electrons are freed, the remaining elec- 
trons experience predominantly the force by the ionic 
background, and they move deep into the potential (see 
their negative values in energy) where they experience 
the full Mie-frequency WMio/t^i = 40/3. In Fig. [THl.d. 
the few electrons with positive energy but small radii 
are those driven back to the cluster by the laser field. 
In Figures [TOfe.f. electrons are strongly aligned (no scat- 
tered points) since the laser field is approaching its max- 
imum (at t = 4-cycles). After the peak of the laser pulse 
fFigs llOb .h') the restoring force of the ions on almost all 
electrons dominates the laser force. 



B. Results for circular polarization 

Since many of the features of energy absorption and 
NLR in a LP laser field are common to the case of CP, 
we only point out the main differences. Equation l|21() 
holds in the CP laser field as well. Figure ITTl shows the 
effective frequency squared vs time for one of the PIC 
electrons, together with the total energy Etot.i{t) and 



2 



solid — (e>^„ltaf 


'' ' / 


dashed — E, , / U 

tot p 


/ 

/" 


/^dashed dot — r / R 


/~ 




/~ 

/ 




' / 




1 ✓ ^ X 

1 







3 3.1 3.2 3.3 3.4 3.5 3.6 

Laser cycles 

FIG. 11: Effective frequency squared (o'eff.i/i^i)'^, excursion 
Ti/R, and the total energy Stot.i = r?(t)/2 — <l>(ri, i) for a PIC 
electron in a CP field vs time in laser cycles. The total energy 
becomes positive only when the NLR is crossed (indicated by 
vertical, dashed line). This result resembles the RSM result 
in Fig. 3] The charge density is p/pc = 40, the peak laser 
intensity is Jo = 2.5 x 10^" W/cm^ 



the excursion / R. One can see that the PIC electron is 
freed (i.e., its total energy becomes positive) only when 
the resonance line (wes/wi)^ = 1 is passed. 

Figure IT^ is the CP analogue of Fig. ^| The results 
are very much similar to the LP case shown in Fig. 1101 
and the arguments made there apply here as well. NLR 
is clearly observed. The main difference is that the PIC 
electrons are nicer aligned towards the resonance point, 
even at early times during the laser pulse (see Fig. 1121b. 
c, d). Almost no scattered particles are visible because 
the dynamics mainly consist of swirling around the clus- 
ter center rather than oscillating through it. The number 
of electrons returning to the cluster is much less so that 
the recombination and rescattering probability is smaller 
in the case of CP. The same is observed in laser-atom in- 
teraction experiments, with important consequences for 
harmonic generation and non-sequential ionization. 

For sufficiently high charge density, the degree of outer 
ionization (defined as the ratio of removed electrons to 
the total number of electrons) is less compared to the 
case of LP laser pulses. The explanation is the same as 
in Sec lII CI the effective laser intensity for CP is by a 
factor of two lower than for LP because of our definition 
of the CP field (p|l. 

One may object that, since the denominator in H21|) 
necessarily increases while the numerator decreases for 
an electron on its way out of the cluster potential, that 
the passage through a point {uul^/uf .Etot/Up) = (2^,0) 
with X some value < (cjMic/'^i)^ is rather the consequence 
of outer ionization than the mechanism behind it. How- 
ever, NLR only occurs at a: = 1, and the results in Fig. 1101 
and Fig.ll2lshow only little spreading along {uics/^\Y 
-Etot = 0. Moreover, the fact that both the single electron 
energies become positive and the radii exceed ~ 2R when 
(woff/wi)^ = 1 indicates that NLR is indeed the responsi- 
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FIG. 12: (color online). Snapshots of PIC electrons in the 
frequency vs energy-plane for CP at times (a) t = 1.5, (b) 
t = 2.0, (c) t = 2.5, (d) t = 3.0, (e) t = 3.5, (f) t = 4.0, 
(g) t = 4.5, and (h) t = 5.0 laser cycles for a laser intensity 
2.5 X 10^"" Wcm"^and (uMic/i^i)^ = 40/3. Other parameters 
as in Fig. |S| The radial positions (in units of R) are color- 
coded. Electrons become free upon crossing the NLR, i.e., 
{Lu!s/uj?,Etot/U^) = {l,0). 



ble mechanism behind outer ionization accompanied by 
efficient absorption of laser energy. 

Finally, Fig. ^| shows the average value of the total 
absorbed energy per electron vs the peak laser intensity 
for cluster charge densities between p/p^ — 3-40 in CP 
and LP laser fields. PIC results are compared with the 
RSM absorption results. The absorbed energy per elec- 
tron in Figs. I13b-d is plotted in units of R'^ujf whereas 
the same results are shown in Figs.ll3b-h in units of the 
ponderomotive energy Up. The PIC results in Figs. 113b- 
d show that the absorbed energy increases linearly in the 
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FIG. 13: (color online). Total absorbed energy per electron 
in units of R^uj^ vs laser intensity for charge densities (a) 
p/pc = 3 (linear resonance), (b) p/pc = 5, (c) p/pc = 20, and 
(d) p/pc = 40. The absorbed energies for LP (red, lighter 
gray) and CP (blue, darker gray) using PIC (dashed) and the 
RSM (solid) are shown. 



log-log representation up to a certain intensity and then 
tends to saturate due to the saturation of outer ioniza- 
tion. One sees that the saturation in the PIC results 
occur close to the RSM threshold intensity. Since with 
increasing charge density the restoring force due to the 
ions increases, the saturation of energy absorption in the 
RSM and PIC occur at higher laser intensities as the den- 
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sity increases from p/pc = 3 to p/ Pc = 40 in Fig. [T^-d. 
When outer ionization and the energy absorption sat- 
urate with increasing peak laser intensity, the average 
absorbed energy per electron divided by U-p (which is 
proportional to the so-called fractional absorption) starts 
decreasing in Figs.ll3b-f. Figures H^.e show that at lin- 
ear resonance p/ Pc = 3, the absorbed energy is already 
high at low values of the laser intensity (< 10^^ Wcm"^), 
and absorption is very efficient as compared to higher 
charge densities p/pc — 5-40, presented in Figs. ICTp-d 
and Figs. I13t-h. In fact, Fig. 113b illustrates that the ab- 
sorbed energy is on the order of ~ lOOf/p (both in the 
RSM and in the PIC) before the saturation of outer ion- 
ization. However, one should bear in mind that at too 
low laser intensities inner ionization would not occur in 
the first place so that in reality there would be no ab- 
sorption at all. 

Energy absorption in CP and LP laser fields at all in- 
tensities and all charge densities are almost equally effi- 
cient. Recalling that with our definition of the CP laser 
field the ponderomotive potential is equal for LP and 
CP while the electric field amplitude is not, we conclude 
that for the absorbed energy Up matters while the NLR 
threshold is determined by the field strength. 

The discrepancy between the PIC and RSM results at 
lower intensities is, again, due to the fact that the RSM 
provides a very well defined potential, although a more 
and more shallow one as the density decreases, in which 
the rigid electron sphere moves. In PIC, the potential 
is initially zero and builds up in the course of outer ion- 
ization. At higher intensities, the potentials for the re- 
maining PIC electrons are closer to the RSM potential 
(see Fig. 01, explaining the improving agreement between 
RSM and PIC results as the laser intensity increases. 

IV. SUMMARY 

In summary, two different approaches to study colli- 
sionless laser energy absorption by clusters, namely (i) 
the rigid sphere model and (ii) particle-in-cell simula- 
tions, were pursued in this work. The goal was to iden- 
tify the dominant mechanism of energy absorption and 
outer ionization of the cluster electrons in near infrared, 



short laser pulses where collisional absorption is known 
to be inefficient. We showed that the cluster electrons 
contributing to efficient absorption and outer ionization 
undergo nonlinear resonance, meaning that the instanta- 
neous frequency of their motion in a time-dependent, an- 
harmonic, effective potential transiently meets the laser 
frequency. Nonlinear resonance is the only possible ab- 
sorption mechanism if the laser pulse is too short for 
the linear resonance to occur (or during the early cluster 
dynamics in longer pulses) and if electron-ion collisions 
(inverse bremsstrahlung) are negligible. In order to prove 
the occurrence of nonlinear resonance we used a method 
to analyze the results obtained from particle-in-cell sim- 
ulations, namely the mapping of the system of electrons 
and ions that interact through their mean ffeld onto a 
system of nonlinear oscillators. 

The occurrence of nonlinear resonance in the particle- 
in-cell simulations presented in this paper resembles the 
nonlinear resonance in the rigid sphere model. For a 
given cluster charge density, there is a threshold inten- 
sity around which the average electron energy displays a 
maximum conversion of laser energy. The threshold in- 
tensity can be calculated using the newly introduced van- 
ishing barrier approximation. The common over-barrier 
approximation — applicable to atoms — fails in the case 
of finite-size potentials and underestimates the required 
laser field strength for ionization. This fact might be 
related to the experimentally observed "enhanced satu- 
ration intensity" in the ionization of finite-size systems 
such as Cgo- 

The efficiency of energy absorption from the laser is 
almost the same for linear and circular polarization, in- 
dicating that "collisions with the cluster boundary" as an 
explanation are misleading. Instead, nonlinear resonance 
is the main absorption mechanism in both cases. To il- 
lustrate this, the rigid sphere model has been extended 
for the circularly polarized laser pulses in this paper. 

Future work will take self-consistent charge state distri- 
butions and mobile ions into account. Preliminary results 
indicate that nonlinear resonance clearly persists under 
these circumstances as well. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft. 
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